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Abstract 



Coupling the galileons to a curved background has been a tradeoff between main- 
taining second order equations of motion, maintaining the galilean shift symmetries, 
and allowing the background metric to be dynamical. We propose a construction which 
can achieve all three for a novel class of galilean invariant models, by coupling a scalar 
with the galilean symmetry to a massive graviton. This generalizes the brane con- 
struction for galileons, by adding to the brane a dynamical metric, (non-universally) 
interacting with the galileon field. Alternatively, it can be thought of as an extension 
of the ghost-free massive gravity, or as a massive graviton-galileon scalar-tensor theory. 
In the decoupling limit of these theories, new kinds of galileon invariant interactions 
arise between the scalar and the longitudinal mode of the graviton. These have higher 
order equations of motion and infinite powers of the field, yet are ghost-free. 
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1 Introduction 



The Galileons [1, 2] are higher derivative scalar field theories with many interesting and 
important properties, including second order equations of motion and novel non-linearly 
realized shift symmetries. Originally formulated in flat space, it is not straightforward to 
couple the galileons to a curved background. Implementing a universal coupling by a naive 
replacement of partial derivatives by covariant derivatives results in a theory with higher 
order equations of motion for the metric. It is possible to add non-minimal couplings to 
restore the second order equations of motion [3, 4], however, doing so inevitably results in 
terms that violate the galileon shift symmetries. 

The brane construction of [5] is an illuminating way to build the galileons 1 . One 
imagines that our space-time is a 3-brane floating in some higher dimensional bulk spacetime, 
and the brane-bending modes then become the galileons. By extending this construction to 
curved bulks and branes, it becomes possible to couple the galileons to different background 
geometries while preserving generalized versions of the galileon shift symmetries, which are 
now associated with isometries of the bulk [7, 8, 9, 10]. The second order equations of motion 
are also preserved in this approach. However, the metric is a fixed background and is not 
dynamical. Making it dynamical corresponds to turning on a zero mode for the bulk metric, 
which breaks the isometries and hence the galileon symmetry [11]. 

On a parallel front, in recent years it has become possible to construct ghost-free (de 
Rham, Gabadadze, Tolley (dRGT)) theories of massive gravity [12, 13] (See [14] for a theory 
review, and [15] for phenomenology review). These theories can be interpreted as the theory 
of a 3-brane embedded in a 3+1 dimensional bulk (i.e. a spacetime- filling embedding), in 
which a dynamical metric is put on the brane [16], and the brane worldvolume action takes 
the form given in [13]. The brane-bending modes become pure gauge modes, or Stuckelberg 
fields, which, in the presence of interaction terms mixing the dynamical bulk metric with 
the induced metric, give the graviton a mass. Interestingly, the Galileon terms emerge in 
ghost-free massive gravity in the decoupling limit [17, 12]. 

In this paper, we combine elements of the brane construction for galileons and massive 
gravity to yield a novel theory that couples a multiplet of scalar fields tt 1 (where J is a 
flavor index, running from 1 to N) to a metric in a way that possesses all three desirable 

1 For a complementary construction of these theories as Wess-Zumino terms, see [6]. 



features: no extra propagating degrees of freedom, a galileon symmetry, and dynamics for the 
metric. We should stress that what we mean by a "galileon" in this work is a generic scalar 
field 7r 7 , nontrivially transforming under the field-space galilean invariance of the theory 



and propagating no extra degrees of freedom than those of a free field. In particular, ir 
does not necessarily have to interact with itself or other fields through the five standard 
galileon terms [2] (or their multi-field generalizations [18, 19, 20]). Moreover, as we show 
below, a certain high-energy ( "decoupling" ) limit of the theory is described by peculiar scalar 
interactions, which significantly differ from the standard galileon interactions, and yet their 
defining properties — ghost freedom and galileon symmetry — are retained. 

In light of findings of Refs. [3, 4], it is not surprising that a theory that can simulta- 
neously achieve these properties is characterized by very special, non-universal couplings of 
the scalars it 1 to the dynamical metric. In particular, as already noted above, the absence of 
higher time derivatives in the equations of motion (or, equivalently, extra propagating ghost 
degrees of freedom) in theories with ir 1 coupled to a massless graviton inevitably results in the 
breaking of galileon symmetries of the flat space theory. Based on the requirements of galilean 
invariance (which arises from nonlinearly realized broken higher-dimensional Poincare sym- 
metries and is therefore automatic in brane constructions), and ghost-freedom, we argue 
below that the metric that can naturally couple to the scalars it 1 describes a massive spin- 
2 field. The essential degrees of freedom on top of the N scalars ir 1 in our construction 
therefore are the five polarizations of a ghost-free massive graviton. 

Formally one can take different points of view towards the theory discussed below. 
On one hand, one can view it as a certain generalization of ghost-free massive gravity, con- 
sistently interacting with a set of scalars ir 1 with the field-space galilean invariance (l.l) 2 . 
Alternatively, one can imagine a four-dimensional effective field theory, obtained by ex- 
tending the brane construction of the galileons to allow for an additional intrinsic metric, 
describing dynamical gravity. Starting with the two non-interacting sectors, consisting of (a) 
the galileons (obtained via invariants formed from the induced metric in the standard way) 
and (b) dynamical gravity, one can ask whether it is possible to construct well-defined mix- 

2 Another generalization of dRGT massive gravity by the "quasi-dilaton" , a scalar realizing a new global 
symmetry in the theory, has recently been considered in [21]. 
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ings/interactions between these sectors, which would not lead to extra propagating degrees 
of freedom. As we argue below, one is quite uniquely led to a generalization of the dRGT 
theories. The byproduct of this construction is that the dynamical metric on the brane will 
inevitably describe a massive graviton, nontrivially transforming under the galileon symme- 
tries. 

Interestingly enough, although governed by similar global symmetries, not all inter- 
actions of it 1 with the longitudinal scalar mode of the graviton in this theory fall into the 
categories of standard galileons or multi-galileons found in [2, 18, 19, 20]. Unlike the usual 
galileon terms and their multi-field generalizations (of which there are a finite number for a 
finite number of fields), these interactions consist of an infinite series of terms with higher 
derivatives, and lead to higher order equations. Yet, these theories are degenerate in the 
higher derivatives, and only two pieces of initial data per field are required to set the dynam- 
ics of the system, so in fact the theory is still free of extra propagating (potentially ghostly) 
degrees of freedom, albeit in a way that differs from the standard galileon interactions (which 
have purely second order equations). 

The paper is organized as follows. We start with reviewing the brane construction of 
galileons in section 2. Sections 3 and 4 describe the ghost-free actions for the induced and 
internal metrics that one can write within this framework, including a generalization of the 
recently proposed zero-derivative interactions, characterizing ghost-free massive gravity. In 
section 5, we introduce the basic model with a flat bulk and comment on its symmetries, 
while sections 6 and 7 deal with the decoupling limit of the simplest version of the theory 
and show how, at least in this limit, non-propagation of extra degrees of freedom is achieved 
despite the presence of higher derivative interactions. Finally, we conclude in section 8. 

2 The general construction 

We begin by presenting the general case of our construction. We will work in arbitrary 
dimension to start, and later specialize to the four dimensional case of interest. The theory 
is that of a (d — l)-brane, with worldvolume coordinates x M , moving in a D- dimensional 
background, D > d, with coordinates X A and a fixed background metric Gab(X). The 
dynamical variables include the brane embedding functions X A (x), D functions of the world- 
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volume coordinates x^. 

We may construct the induced metric g^(x) via 

r) X A c)X B 

9M = ^r-^r G ^(x(x)) . (2.1) 

In addition to the induced metric, there are other geometric quantities associated with the 
embedding, such as the extrinsic curvatures and the twist connection (see Appendix A of [20] 
for a complete description of these quantities). 

We would like the action on the world- volume to be invariant under re-parametrizations 
of the brane x^ — > x^ — ^(x), under which the embedding functions are scalars, 

5 g X A = ed,X A . (2.2) 

The induced metric (2.1) transforms as a tensor under these gauge transformations, 

Sgg^u = C^g^ = ^dxg^v + g Xv + d v £ x g^ x ■ (2.3) 

Gauge invariance requires that the action be written as a diffeomorphism scalar, F, of the 
induced metric g^, its covariant derivatives V M , its curvature R p c , and the other induced 
quantities such as intrinsic curvature and twist (which we denote with ellipses), 

S = J d d x y^F (g^, V„, R p a , u , ■■■) . (2.4) 

In addition to the gauge symmetry of re-parametrization invariance, there can also be 
global symmetries. If the bulk metric has a Killing vector K A (X), satisfying the Killing 
equation 

K c d c G AB + d A K c G CB + d B K c G AC = , (2.5) 
then the action will have a global symmetry under which the embedding scalars X A shift, 

6 K X A = K A (X) . (2.6) 

The induced metric (2.1) and other induced quantities, and therefore the general action (2.4), 
are invariant under (2.6). 

We may completely fix the re-parametrization freedom (2.2) by fixing the unitary gauge 

= xf, X\x) = tx\x) . (2.7) 



In this gauge, the world-volume coordinates of the brane are identified with the first d of the 
bulk coordinates. The remaining unfixed fields, vr 7 (x), I = 1 ■ ■ ■ N, where N = D — d is the 
co-dimension of the brane, measure the transverse position of the brane. The gauge fixed 
action is an action solely for ir, 

Sg = J d d x V^F (g, v ,\7„R a ^ u ,- ■ - )\ x ^ XI=nI . (2.8) 

The form of the global symmetries (2.6) is altered once the gauge is fixed, because the 
gauge choice (2.7) is not generally preserved by the global symmetry. The change induced 
by K A is Skx^ 1 = K^(x,ti), Sr^ 1 = K\x, it), so to maintain the gauge (2.7), we must 
simultaneously perform a compensating gauge transformation with the gauge parameter 

C mp = -i^(*,7r) • (2-9) 
The combined symmetry acting on the fields it 1 is now 

(5 K + 5comp)7r 7 = -K»(x, n)d^ + K\x, tt) , (2.10) 
and is a global symmetry of the gauge fixed action (2.8). 

In addition to the induced metric (2.1), we now introduce an additional world-volume 
metric g^ v (x) onto the brane. We demand that this obey the same transformation laws as 
g^v, and so we declare that it is invariant under the global symmetries (2.6), but transforms 
clS db tensor under (2.2), 

5 K g»u = 0, (2.11) 
5 g g^ u = £>i9nv = £ X d x g^ + d^ x g Xv + d^ x g^ ■ (2. 12) 

We are now free to add to the action terms which are diffeomorphism scalars constructed 
from the intrinsic metric g^ v and its associated covariant derivative and curvature, 

S g = J d d x 4^gF 0^, V M , i?V) > ( 2 - 13 ) 
as well as terms that mix the intrinsic metric with the induced metric and other quantities 3 , 
S mix = J d d x ^~gF (g^ V M , R a ^ u , g, u , V„, R a ^ v , ■ ■ ■ ) . (2.14) 

3 We have chosen \J—g to be the universal measure factor for each term in this part of the action. This 
entails no loss of generality, since this choice may be traded for y/—g, or even a geometric mean such as 
(— g) 1 / 4 (— g) 1 / 4 , by pulling out or absorbing powers of factors such as det (<7 _1 g) into F. However, whatever 
choice is made the absence of the ghost should be investigated, as we will do in subsequent sections. 
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Once the unitary gauge (2.7) is fixed, the fundamental fields of the theory are g^ v and 
tt 1 , and the global symmetries act as 

8 9tlv = -K x (x, Tt)d x9ilv - d, L [K x (x, tt)] gxu - d u [K x (x, tt)] 9tlX , (2.15) 
H 1 = -K*{x, 71)8^ + K\x, tt) . (2.16) 

(Here, we must act by the derivatives on the argument of ir(x) within K^, i.e. [K x (x, 7r)] = 
d^K x + t^p-d^ir 1 .) Note that in the unitary gauge the intrinsic metric g^ v transforms non- 
trivially under the global symmetry, due to the compensating gauge transformation (2.9). 
Under this (in general non-linear) transformation, the induced metric (2.1) transforms in 
precisely the same way as the intrinsic metric (2.15), 

= -K x (x, n)d x g^ - d„ [K x (x, tt)] g Xv - d v [K x (x, tt)] g„ x . (2.17) 

This construction allows us to have scalars with galileon-like shift symmetries given 
by (2.16), coupled to a dynamical metric, which now carries a non-trivial transforma- 
tion (2.15) under the galilean symmetries. It remains to ensure the final desired property - 
that the action is free of ghosts. 

3 Ghost-free actions 

If, as for the original galileons, the actions are to be free from extra Boulware-Deser-like [22] 
degrees of freedom, they must take a specific form. For the part of the action S g , depending 
only on the dynamical metric g^, we know that the only possibilities giving second order 
equations of motion are the Einstein-Hilbert term, the cosmological constant, and the higher 
Lovelock invariants [23, 24] if the brane has dimension d > 4, 

S 9 = ^ J d d xV^g[-2A + R[g} + ---} • (3.1) 

For the term Sg, depending only on the induced metric g^, the possibilities are the 
Lovelock terms of g, as well as the boundary terms associated with Lovelock terms in the 
bulk, as detailed in [20]. For co-dimension 1, these are the Myers boundary terms [25]. In 
the unitary gauge (2.7), these lead to the galileon terms for the single tt field [5]. For higher 
co-dimension, the surface terms are more limited and difficult to catalogue [26, 27]. In the 
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unitary gauge they lead to the multi- field galileons [20]. In all cases, the leading term is the 
DBI term for the induced metric, which contains the kinetic term for the n 1 fields, and so 
we write this part of the action as 

S s = -T J d d x V=g~+--- , (3.2) 

where T is a constant of mass dimension d, and the ellipses denote the possible higher order 
Lovelock and boundary terms. 

For the mixed terms, it is not immediately obvious what the most general ghost-free 
terms are. However, if we restrict to terms depending only on g^ y and g^, with no higher 
derivatives, we can take a clue from the dRGT theory [12, 13] of massive gravity and the 
related models of bi-gravity [28], all of which have been shown to be ghost free [29, 30, 31, 
32, 33, 34, 35]. These models contain interaction terms between two metrics — the second 
metric is a fixed fiducial metric in the case of massive gravity and a dynamical second metric 
in the case of bi-gravity. In this paper we will choose the form of the interactions to be the 
same, but the second metric will be the induced metric (2.1), containing the ir 1 degrees of 
freedom. The interactions can be constructed through the tensor K) l v = 8** u — \/ g^gxu, in 
terms of which the relevant piece of the action is given as follows 

Smix = -"Tp / d^x^g ^^2a n S n (JC) , (3.3) 

J n=2 

where S n (M), < n < d, for a d x d matrix M^, are the elementary symmetric polynomials 4 

S n {M) = M^\...M»;l (3.4) 

and y 1 g~ x g is the matrix square root of the matrix g hLa g al/ . The 0^4 are free coefficients, 
while «2 = —8 is required for the correct normalization of the graviton mass. It is often 
convenient to work in terms of the expanded action 

yri r m 2 d 

s mi , = / d d xV^ T J2^s n (Vg = ^) , (3-5) 

J n=0 

where /3 n can be expressed in terms of the two free parameters 0:3,4. We will use the latter 
representation of the mixing terms below. Note that the n = and n = d terms in the 

4 Our anti-symmetrization weight is [/ii . . . fi n ] = ^t(mi ' ' " A*n + • • • )• The appearance of the symmetric 
polynomials and their relation to the absence of ghosts can be naturally seen in the vielbein formulation of 
the theory [35]. See appendix A of [35] for more details on the symmetric polynomials. 
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latter sum are redundant, since these reproduce the cosmological constant y/—g and DBI 
term \f—g respectively. 



4 Co-dimension zero: massive gravity 

Our construction contains ghost-free dRGT massive gravity as a special case. When the 
co-dimension is zero, we are embedding a (i-dimensional world-volume into a bulk space of 
the same dimension. The fixed bulk metric G^ u therefore has the same dimension as the 
brane metric. 

In unitary gauge there are no 7r fields, and the induced metric is the bulk metric, 

9^{x) = G^(x) , (4.1) 

so that the global symmetries are the Killing vectors of G^ v (x): (, x d\G^ v + d^ x G\ u + 

d^G^x = 0. The intrinsic metric then transforms linearly as a tensor, 

Sg„u = i X d x g„u + d^ x g Xu + 8^ x g„ x . (4.2) 



The action Sg contains no dynamical variables and can be dropped. If the bulk metric 
is flat, G^ = rin U , the action S g + S m ^ is precisely the Lorentz invariant dRGT massive 
gravity of [12, 13], in the unitary gauge. Further, the Poincare invariance of these actions 
comes from (4.2). For a general bulk metric, the theory is that of massive gravity with a 
general reference metric [30], and the global symmetries (4.2) are precisely the isometries 
of the reference metric. These theories are all ghost free [29, 30, 31], meaning that they 
propagate, non-linearly, precisely the number of degrees of freedom of a massive graviton 
and no more. 

Away from the unitary gauge, we have 

f)vp f)~x a 

which is nothing but the Stiickelberg replacement used to restore diffeomorphism invariance 
to massive gravity [36, 16]. 
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5 Flat bulk case 



We now return to general co-dimension N, but specialize to a flat bulk metric Gab — Vab- 
The isometries are the Poincare transformations of the bulk, 

5 P X A = K A (X) = u A B X B + e A , (5.1) 

where e A and the antisymmetric matrix u A B are the infinitesimal parameters of the bulk 
translations and Lorentz transformations respectively. 

The unitary gauge (2.7) is not in general preserved by the Poincare transformations, but 
the gauge is restored by making the compensating gauge transformation, SgX^ = t^dyX^ = 
with the choice 

Cmp = - «V - C . (5.2) 

The combined transformation 5 pi = 5p + 5 g then leaves the gauge fixing intact and is a 
symmetry of the gauge fixed action. This symmetry acts on the remaining fields as 

Sp^ 1 = -ou^Xdy - e^ir 1 + + uif^d^ 1 + e 1 + u^tt j , (5.3) 

where the first two terms in this expression are unbroken spacetime rotations and transla- 
tions, respectively. The second two terms are a DBI symmetry corresponding to the broken 
boosts in the extra dimensional directions (which becomes the galileon symmetry for small 
ir 1 ). The fifth term is a shift symmetry, corresponding to broken translations into the trans- 
verse directions. Finally, the last term is the unbroken SO(N) symmetry in the transverse 
directions, which appears as an internal rotation among the ir fields. The symmetry breaking 
pattern is ISO(l, D - 1) -> ISO(l, d - 1) x SO(N). 

The induced metric in unitary gauge is 

g^v = rjpv + d^dv-Ki , (5.4) 

and using (2.15), we can determine how the global symmetries extend to the metric. The 
metric transforms linearly, as a tensor, under the unbroken d- dimensional Poincare symmetry, 
and is invariant under the unbroken SO(N) internal symmetry. The broken DBI shift 
symmetries, on the other hand, extend non-trivially to the dynamical metric, 

&9iu> = ^/ A?r/ d\9ii» + ufdpir 1 g Xv + w/^tt 7 g^ x . (5.5) 
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This is a diffeomorphism with parameter £ M = uj^tt 1 , since its origin was nothing but a 
compensating diffeomorphism to restore unitary gauge. The induced metric (5.4) transforms 
in the same way, as 

ftfjfiv = ufit 1 dxg^u + ^/<9 /t 7r 7 g\u + u^d^n 1 . (5.6) 

In summary, in the unitary gauge the galileons 7r 7 and the intrinsic metric transform 
as tensors under the unbroken d dimensional Poincare symmetry, as a vector and singlet 
respectively under SO(N), and as follows under the broken galileon symmetries, 

Six 1 = u 1 ^ + uf^dy + e 1 , 

5g^ w = ufir 1 d x g^u + co^d^ir 1 g\ v + uj^dyTi 1 g^\ . (5.7) 

This is how the DBI galileon symmetry extends to the metric. It is a global symmetry. 
In the unitary gauge we are working in, there is no diffeomorphism invariance (assuming 
there are interaction terms S m ix in the action). As a consequence, the graviton described 
by g^ u will be massive. In this sense, it is natural for the galileons to couple to a massive 
graviton. 

If we choose not to go to unitary gauge, diffeomorphism invariance on the brane then 
remains intact. In this case we have the Stuckelberg fields X^, and the induced metric takes 
the form 

dXP8X° a Ia 
9fiu = dx» dx" Vpa + At?r u1Tl ' ( ) 

Because we now still have diffeomorphism invariance, the induced and intrinsic metrics are 
invariant under the global symmetries (5.1), and transform as tensors under diffeomorphisms. 
The fields X^ and it 1 are scalars under the diffeomorphisms, and transform together as (5.1) 
under the global symmetries. 



6 Small field expansions and decoupling limits 

If the action has a background solution g^ = rj^ for the intrinsic metric, we may expand 
about it in fluctuations h^, 

g^v = + . (6.1) 
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The unitary gauge non-linear transformation laws (5.7) expanded around this background 
are then 

5tt j = u 1 ^ + ufir J d IJi ir I + e T , 

Sh^ = ui^du-n 1 + uivd^ 1 + w/V d\h^ u + u^d^ir 1 h Xu + u^d^n 1 h^x , (6.2) 

and we see that the metric fluctuations must transform along with the galileon fields. 
Defining the fluctuation around the induced metric via 

Hf.u = g^u - 9nv , (6.3) 

we have in unitary gauge 

= /v - di^ l d v T\i . (6.4) 

To perform the Stiickelberg expansion, we simply leave the gauge unfixed, so that the 
induced metric takes the form (5.8). The fluctuation (6.3) can then be written as 

H^ u = + 7]^ - d^dyX^r}^ - d^d^i . (6.5) 

As in massive gravity, we can then introduce another Stiickelberg field <ft to deal with the 
longitudinal mode through the following replacement, by expanding X p around its unitary 
gauge value, 

X p = x p + A p - rTd a <t> ■ (6.6) 
The action then has the infinitesimal gauge transformations, 

Shfj, v = <9 M £„ + d u ^ + C^h^ , 
5(f) = A , 

Sn 1 = Cd^ 1 . (6.7) 

As is usually done in massive gravity, we ignore the vector mode A 11 which carries the helicity 
one components of the massive graviton at high energy (a consistent truncation), since these 
do not generally couple to matter at the linearized level. Putting them to zero is consistent 
with the equations of motion. Moreover, there is an enhanced U(l) symmetry for this vector 
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in the decoupling limit that guaranties that it propagates two degrees of freedom. We then 
have the replacement 



H^ u = hfa, + 2d^d„(p - d^d x <fi d u d X (p - d^d^i . (6.8) 

For generic choices of the action, these theories describe a massive graviton coupled to 
the galileon fields tt 1 , with coupling such that the unitary gauge action is invariant under the 
galileon symmetries (5.7). Away from the unitary gauge on the other hand, the longitudinal 
mode of the massive graviton is described by the scalar 0, which appears in addition to the 
n 1 when we restore the diffeomorphism invariance by not fixing unitary gauge. 

For definiteness, we now focus on d = 4. The action is 



/M m f > 

d 4 xl -j 1 ^ R[g]-—Y,PnS n (yg^Zj^ 

I L n=0 J j 

(6.9) 

Here Ao and f3 n are order one dimensionless constants, independent of the mass scales m and 
M-p\. (We have chosen the mass scalings of the various terms so that there will be an inter- 
esting decoupling limit, with new ingredients beyond those appearing in the corresponding 
limit of massive gravity. ) The ellipses in the final term denote the possible higher-order Love- 
lock and boundary terms, each of which has its own independent order one coefficient and 
is suppressed by the mass scale m (which will ensure that non-trivial galileon interactions 
survive the decoupling limit). 

The coefficient 04 is redundant with Ao, and so we set 04 = 0, and expand around flat 
space in the unitary gauge g^ u = rj^ + h^, g^ v = r]^ v + d '^ii 1 d u ir j . Tadpole cancellation, 
ensuring that flat space is a solution, requires 

A, + 30i + 3f3 2 + p 3 = 0. (6.10) 

At quadratic order, we find a Fierz-Pauli massive graviton. One of the /3's can be absorbed 
into m 2 , which we do by demanding 

8 -ft -2ft -A, = 0. (6.11) 

This ensures that the Fierz-Pauli massive graviton with a mass m propagates at quadratic 
order. 
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In addition we find a kinetic term for the it fields, 

J\ 2 



2 P1 



Ao + I (A + 3/3 2 + 3/3 3 ) 



(dyr 7 ) . (6.12) 



Thus, provided A + | (/Si + 3/9 2 + 3/9 3 ) > 0, the theory propagates, in addition to the massive 
graviton, iV healthy (having the correct sign kinetic term) scalars on this background. 

In the end, we have four free parameters for the flat space theory (plus those correspond- 
ing to the higher Lovelock terms): the graviton mass, two independent /3's corresponding 
to the two free parameters in the interactions of dRGT massive gravity, and a remaining 
independent parameter which corresponds to the strength of the 7r kinetic term. Note that, 
if we had been looking for curved (A)dS solutions, there would have been an additional 
parameter corresponding to the curvature of the background. 

We now make the replacement (6.8) and expand in powers of the fields. After canonical 
normalization of the various kinetic terms via 

h ~ Mp\h , <p ~ m 2 Mp\(j) , it ~ mMpiir , (6.13) 

one can examine the interaction terms to determine their associated interaction scales. 

First focus on those interaction terms arising from S m i x (or from the DBI term in Sg). 
By virtue of the Stiickelberg replacement (6.8), <fi always appears with two derivatives, it 
appears with one derivative, and h appears with no derivatives. A generic term with rth 
powers of h^ u , powers of n 1 and powers of <fi, reads 

~ m 2 M^h nh (dTr) n -(d 2 (f)) n ^ ~ A'" nh " 2 ^" 3n ^ nh (97r)^(9 2 0)^ , (6.14) 

where the scale suppressing the term is written as 

A A = (M P1 m ) , A = ■ ■ — . (6.15) 

+ n n + n h - 2 

Since we always assume m < Mp\, the larger A, the smaller is this scale. Note that 
must be even, by virtue of the way it enters in (6.8), and we have + ua + > 3, since 
we are only considering interaction terms. The terms suppressed by the smallest scale are 
self-interaction terms, = rih = 0, which are suppressed by scales > A 5 and < A 3 . 
These terms, however, all cancel up to a total derivative due to the special structure of the 
ghost-free massive gravity interactions [12, 13]. 
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The scale A3 = (Mpjm 2 ) becomes the lowest scale, and is carried by terms of 
schematic form {n > 1) 

~ TZ=lk&4>) n > (6-16) 

A 3 

and 5 

~ ^f^T • (6-17) 

(Note that here we have included the term mixing h and for 71 = 1, even though it is 
a kinetic mixing term and not an interaction term, because it is from this mixing that <ft 
acquires its kinetic term and its canonical normalization.) All other terms carry scales higher 
than A 3 . There are a finite number of terms of the first type (6.16), and they take the same 
form as they do in massive gravity [12, 13]. However, there are an infinite number of terms 
of the second type (6.17). 

We now return to the possibility of higher Lovelock terms in Sg. In the unitary gauge, 
before any decoupling limit, these are the DBI galileons [5, 37] for N = 1, and their multi- 
field generalizations for higher N [20]. As we have mentioned, these terms are suppressed by 
the scale m. For example, for N = 1 the leading Lovelock term beyond the DBI kinetic term 
is the trace of the extrinsic curvature, which leads to a cubic DBI galileon in the unitary 
gauge action, 

S- g = -M^m 2 f d A x (\ 0V fTTW¥- Irr^d^d^d-Ti + •••). (6.18) 

J V m 1 + [airy J 

Away from the unitary gauge on the other hand, restoring the Stiickelberg field and canon- 
ically normalizing via fr ~ mMp\7r, the various galileon terms yield interactions of the form 

2 2 fd 2 4>Y' ( '0* Y" fd 2 *\" , s 

~ M - m \i) \m^) (at) • »-- 1 >» = ». 1 .-,4. (6.19) 

The terms with multiple powers of the the factor ^ arise from expanding out the square 
roots and denominators of (6.18). The interactions suppressed by A3 are those with n n = 
n + 2. Of these, the ones with n§ = are precisely the terms which survive the limit which 
recovers the normal galileons from the DBI galileons [5]. Those with > 0, of which there 
are an infinite number, describe the coupling of the galileons with the longitudinal mode of 

5 Had we not neglected the vector mode A^, we would have seen that terms of the form dAdA(ddir) n /A 3n 
also survive in the limit at hand. 
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the graviton. We work out these couplings for the case of the cubic galileon in Appendix A. 
All other terms are suppressed by scales larger than A 3 . 

If we take the decoupling limit, 

Mpi -> oo , m -> , A 3 = (M p m 2 ) 1/3 fixed , (6.20) 

then all interactions with scales greater than A 3 are set to zero, and the DBI-galileon terms 
become the normal galileons. The only part of the Einstein-Hilbert action that survives the 
decoupling limit is the quadratic part, 

r2 



/M l 
d*x -^h^e^haf, , (6.21; 



where the kinetic operator for the graviton is that of linearized Einstein gravity 6 . The gauge 
symmetries (6.7) in the decoupling limit become the linearized versions of those considered 
above 

Sh^ = d^y + , 5A IM = d^A , 
5<p = 0, 5vr 7 = 0. (6.23) 

As an explicit example, consider the choice of coefficients f3 = —24, /?i = 8, f3 2 = /3 3 = 
A = so that the DBI term is gone, and set to zero all the higher Lovelock terms within Sg 7 . 
This choice corresponds in pure massive gravity to the model which contains no nonlinear 
scalar-tensor interactions in the decoupling limit [12]. The corresponding dRGT action, 



written in terms of the quadratic, cubic, and quartic terms in the matrix /C^ = b yi v — \J g^ x g\ v 
and its traces [13], can be rewritten for these particular coefficients in the form given in [38] 

/M 2 r / \ i 

d 4 x^ v ^[ J RM+2m 2 (3-Tr v ^)J , (6.24) 

which is often referred as the "minimal model". Unlike in pure massive gravity though, 
here g is going to be an induced metric on the brane. Thus the interactions we find will be 
entirely due to the extension of the theory we have developed here. 



Explicitly, 

S^p = (firfl - rTVap) □ - 29^8^ + d»d»r) a(S + d a d p rf '. (6.22) 
7 This corresponds to the choice C3 = ds — — in the notations of [12, 14]. 
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The decoupling limit Lagrangian is 8 



C 



M 



dec 



a8 



'PI' 



V (rr - 



1 



2 <?£ - d»d v i 



(6.26) 

Diagonalizing via /i^ = fa' + m 2 077 MJ/ , we find a free decoupled graviton, and coupled 
interacting scalars, 



C 



Mi 



dec 



— -m 2 (<90) 2 — -—r, - 

4 v V) 2 5i-d» t d v (, 



(6.27) 



The terms involving n 1 are new to this model, and do not appear in pure massive gravity. 
Note that, unlike massive gravity, there are an infinite number of scalar interaction terms 
that survive the decoupling limit. 



7 Equations of motion and ghosts 

The equations of motion obtained from (6.27) are not second order. To see this, we need 
only expand to cubic order in the fields, £ cu bic ~ d fJj d u (f>d tJ "rT I d u 7i I . The equation of motion, 
for example, is third order, ~ d^d v [d^n 1 d u ni) . 

Higher order equations are generally associated with extra ghostly degrees of freedom. 
In dRGT massive gravity, the decoupling limit is second order and contains no extra degrees 
of freedom, as it must since the entire theory has no such extra degrees of freedom. The 
higher order equations we are finding here are naively worrisome, because if the decoupling 
limit contains extra degrees of freedom, the entire model is not ghost free. 

However, higher order equations do not necessarily imply the existence of extra degrees 
of freedom. As we will now show, the Lagrangian (6.27) in fact contains no additional ghostly 
degrees of freedom, despite the higher derivatives. 

8 For deriving the decoupling limit, it is convenient to write the Lagrangian in terms of the tensor K.^ v = 
K - vV A Sa^ C = ^V^9 (■%] + 2 ™ 2 (-1 + £^)) • Putting in the Stiickelberg helds, 

KP V = - ^&%-g^H Xv , = V + 2IV - K» - P ^ > ( 6 - 25 ) 

with Il^y = d^d u (f), and P^ u = d^ir 1 d^nj . We then use the relation (in matrix notation with ( ) the trace) 
(SIC) = — — II) _1 P) to find the terms involving tt, and JC = 11+ \h— jhU — jHh for the terms involving 
h. 
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7.1 A toy example 



As a warmup to proving this, consider first the simpler + 1 dimensional version of the scalar 
part of the action (6.27), 



,5' 



1 .. 



This is a higher-derivative Lagrangian, and the equations of motion are naively fourth order 
(third order if we expand to cubic order in the fields, as we did in the full case above), 



5S 

ss_ 

~54> 



d 
dt 



7T 



+ 2 



7T 



(7.2) 
(7.3) 



.1 + 0) 2 . 

As with the full theory, this seems worrisome, since it raises the possibility of extra degrees 
of freedom which are ghosts. 

However, the number of initial data needed to solve this system is only four, consistent 
with there being only two degrees of freedom. To see this, note that the ir equation implies 
is a constant, which when substituted into the equation implies = 0. Substituting 
this back into the tt equation then yields if = 0. Thus the equations above are, in fact, 
equivalent to the free field equations = if = 0, and there are therefore no extra degrees of 
freedom. 

Note that Ostragradskii's theorem [39] does not apply to (7.1) (or to the full model (6.27)), 
since one of the conditions of the theorem, that the Lagrangian be non-degenerate in the 
higher derivatives (i.e., the matrix obtained by variation of the action with respect to second 
derivatives be nondegenerate) , is not satisfied. 

The absence of extra degrees of freedom can also be seen directly at the level of the 
action. Starting with (7.1), we introduce an auxiliary field a to render the action polynomial 
in the fields, 



S 



dt 



1$> _ + 

T 2 v 



7TO" 



(7.4) 



which when substituted back into (7.4) 



The equation of motion for a is then a 
recovers (7.1). Integrating by parts to remove the second derivatives from 
equivalent first order action, 



,5' 



dt 



1 12 1 2 

-0 - -a +aa(, 



no 



we find an 



(7.5) 
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We now Legendre transform to find an equivalent Hamiltonian action. The canonical 
momenta are 

p^ = 4> + a&, p n = 0, p a = a<f> - 7r . (7.6) 
There is a primary constraint 

P, = , (7.7) 
for which we introduce the multiplier A, and the action then takes the form 

>2 



S 



dt 



+ p w TT + p a <J 



1 , » 2 Pd> , \ a \ 

(p° + n ) + — (p° + 71 > + y ) ~ 



(7.* 



It is now straightforward to see that the primary constraint (7.7) generates a secondary 
constraint, leaving four phase space degrees of freedom, or two Lagrangian degrees of freedom. 

Alternatively, we may solve the primary constraint directly in the action, 



S 



dt 



+ Pa<? 



1 



2a' 



(p a + 7T y + ^( P(7 + 7r) + "- 



a 



Now it is an auxiliary field and can be eliminated through its equation of motion, n 
leaving 



S 



dt 



+ Pa<j 



1 2,1 2 

2 p + + r 



(7.9) 

IP-Pa, 
(7.10) 



Thus, we see explicitly that there are exactly two degrees of freedom, since renaming p a — > q, 
a — > —p q , the action is equivalent to that of two free particles with positive energy, 



5 



dt 



+ PqQ 



± 2 , x 2 

2 P + + 2 



(7.n; 



7.2 No extra degrees of freedom in the decoupling limit 



We now apply apply a Hamiltonian analysis to the full scalar action in the decoupling limit 
of the theory (6.27). Setting 3m 2 /4 = 1/2, and specializing to the case of a single n field for 
simplicity, the action (6.27) we are studying is proportional to 

1 ,„ ,,9 1„ 1 

2 



,5' 



d A x 



1 

" 2 Ufl " 5t - d»d v 



-d V TX 



(7.12) 



To eliminate the inverse powers of derivatives and work with a local action, we introduce 
an auxiliary vector field and write the following equivalent action, 



,5' 



d 4 x 



(7.13) 
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Integrating out f2 M through its equations of motion recovers (7.12). Now we make a (3 + 1) 
decomposition of the Lorentz indices, and do some integrations by parts to remove all the 
double time derivatives from d>, 



S 



d A x 



l X2 1 



2W) 2 -\{^? + \m 2 +n°n°< 



d x 



(7.14) 



As in the toy model, we now have a Lagrangian which has at most first time derivatives, 
so we may pass to a Hamiltonian form of the action in standard fashion. The conjugate 
momenta are 



P<f> 



n°<p + 7T , 
, 
. 



From this, we see that we have the primary constraints 

Pit = , PQi = . 

The Hamiltonian density, on the constraint surface, is 
% = P<t><fi + Pno&° +PkK + Po.i& - C 



(7.15) 
(7.16) 
(7.17) 
(7.18) 



(7.19) 



;(PfP -tt) - 



(7.20) 



Thus, solving the constraint inside the action gives the following equivalent Hamiltonian 
form of the action 



s 



d x 



(7.21) 



We can now see that the system describes precisely two fields' worth of degrees of freedom: 
both 7r and fl l appear algebraically and can be integrated out, leaving a Hamiltonian action 
depending only on the phase space variables p^, 4>,Pqo, Q°, for which there are standard 
unconstrained first order Hamiltonian equations. 
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Doing this explicitly, we first integrate out n with its equation of motion 



7r = Pn o - ft% + ft u ft l <9ift u , (7.22) 
which, substituted back into the Hamiltonian density, yields 

+ i(ft ) 2 - i(ff) 2 + l^iPdidjff) . (7.23) 
Next we eliminate ft* through its equation of motion 

ft* = (A -1 ) iJ '(%>n° - PA^°) , (7-24) 
where A*- = 6 l , — (d l O,°)djO,° — d l dj(p , giving the Hamiltonian 

n = y + -Mft )^- 1 )^^^ -p^-ft ) + ^(V0) 2 + ^(ft ) 2 . (7.25) 

This describes two fields with nonlinear, spatially non-local interactions. 
At the linear level, we have 



U = f + ^(Vpno) 2 + i (V0) 2 + l(ft ) 2 , (7.26) 

so if we redefine p n o = x an d ft = — p x , we can see that this describes two ghost-free, 
non-tachyonic modes 

2 2 

H = ^ + ^ + i(V0) 2 + i(Vx) 2 . (7.27) 

As already noted above, we have assumed that the vector mode also present in the 
decoupling limit does not get excited 9 , so that we may set = 0. This is a consistent 
truncation of the action, since only enters at the quadratic level in the action. (It is 
also consistent quantum mechanically, in the decoupling limit, due to the Z 2 symmetry 
A^ — > — A^.) For arbitrary excitations of this mode, our Hamiltonian treatment has to be 



9 Note that the presence of the infinite number of interactions of the form dAdA(ddTr) n in the decoupling 
limit is also essential for the (nonlinearly realized) invariance under the broken bulk Lorentz generators in 
(5.1). While the inhomogcneous piece in the galilcon transformation is 6ir = us^x 11 , the vector Stiickclbcrg 
mode shifts under this generator as SA^ = —ui^n. The infinite number of these terms then should relate by 
this symmetry to the infinite number of scalar interactions found in (6.27). 
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modified. We have no handle on the infinite number of dAdA(ddir) n interactions present in 
the decoupling limit (even in dRGT gravity their form is in general not known beyond the 
cubic order [17], see however [34, 40]). However, it is plausible to expect that the presence of 
the vector mode does not spoil the ghost-free property of the decoupling limit because of the 
enhanced U(l) symmetry of the limiting action; this is precisely what happens in ghost-free 
massive gravity. 

If curvature invariants composed of the induced metric are added as implied by an 
ellipsis in (3.1), the proof of ghost freedom becomes increasingly difficult; the special structure 
of the resultant decoupling limit interactions however leads us to conjecture that the absence 
of the extra degrees of freedom carries through to this well. The simplest such 

interaction, corresponding to a cubic galileon for 7r in the decoupling limit, is considered in 
Appendix A. 

8 Summary and prospects 

We have introduced a model which couples a scalar ir 1 to a dynamical metric in a manner 
which respects the galileon symmetries. The metric to which the galileons couple is a massive 
graviton. The model can be considered as an extension to higher co-dimension of ghost-free 
dRGT massive gravity, or as an extension of the brane construction of the galileons where a 
dynamical metric is added to the brane. 

We have derived explicitly the decoupling limit of the model around flat space, for a 
specific choice of parameters, and have shown that there are no ghosts. We have not proven 
that there are no ghosts beyond the decoupling limit, though we expect that there should 
not be, since the model is based on the ghost-free constructions of dRGT massive gravity 
and galileon theories. 

This model should provide a completely consistent framework within which to investi- 
gate the implications of galileon invariance in, for example, cosmology. 
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A Adding a cubic galileon 

In this appendix, specializing to the case of a single extra dimension for simplicity, we 
derive the decoupling limit interactions resulting from adding the extrinsic curvature term 
Ml x m j d 4 Xy/=g K(g) to the r.h.s. of (3.1). 

As above, we we work in the gauge = X 5 = n(x), so that the function 

= 7r (x(X)) — X 5 = defines the embedding. The vector ua, normal to the brane has 
the following components, 

9a® d,n 1 , An 

UA \r) AB d A <S> d B <S>\V 2 M (l + a-7r<9 a 7r)V2' U5 ~ (1 + ^tt^tt) 1 /* ' 1 ' 

where the operator 8 denotes differentiation with respect to the bulk coordinate, 

d u = ^— -— = (d X^) = Af(x)d p . (A.2) 

dX" Oxp v p 1 dxP y ! 

The trace of the extrinsic curvature is then given by 



A R AR 1 ( ',=, 8 i/K 8 V TT 8^ 8 V Tl \ 

K = 4 <S r = n A , B n» = + (□* - TTs^- ) ■ (A ' 3) 



The last step is to evaluate \/—g, 

= d^X a d u X^ aP + d^d v ix , (A.4) 
which, multiplied by two factors of the operator A on both sides, gives 

det {A^A p % v ) = det ( VXp + 8 x n8 p n) = -1 - (8tt) 2 g = - (l + (<9tt) 2 ) det (A~ 2 ) , 
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so that we have 



10 



M&m I d\ V^K = M^m I d 4 x det(dX) ( Chr - ^if^l ) . ( A . 5 ) 



1 + d a nd a -K 

In the weak field, and the decoupling limits, this leads to an extra cubic galileon with the 
5D derivatives d — > 8 instead of ordinary ones in (6.27), 

£ dcc D -Ml x m J d 4 x det (1 - dd<p) d^d^Wit , (A.6) 

where (in the decoupling limit) <9 M = [ (1 — dd^)" 1 ■ d\ . In (0 + 1) one dimension, the latter 
expression becomes a surface term, as can be seen by reparametrizing the time coordinate 
t-tt = J dt(l-<j>). 
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